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1. Introduction 



The study of classical kinetic models in one dimension has revealed interesting physical 
properties (e.g. non-equilibrium phase transitions) and mathematical structures. Moreover, 
one finds close connections with quantum-mechanical spin problems. A prominent example 
is the diffusion of particles with hard-core repulsion on the sites of a chain which is coupled 
to reservoirs at both ends [1 — 4] . This model shows at least three phases which differ in 
their density profiles and the current through the system. Boundary effects, i.e. the rates 
at which particles enter and leave the system, play an essential role. Mathematically, the 
model can be described as a spin one-half problem and, if the dynamics consists of single 
particle processes, the time evolution operator of the master equation has the form of the 
Hamiltonian of the Heisenberg model with boundary fields at both ends. The stationary 
state can be written in the form of a matrix product, where the weight of a configuration 
is given by an expression of the form {W\ABBAB ■ ■ .\V) with operators A, B and vectors 
\V) in an auxiliary space. Such states, which generalize simple product states so 
as to give non-trivial correlations, were first found in a problem of lattice animals [5] 
and for certain quantum spin chains [6, 7] . They have also been encountered in diffusion- 
coagulation models [8] . The detailed mechanism is somewhat different in each of the three 
cases. 

In the following we will study this model with a more general type of dynamics. The 
time is taken discrete and in each time step, hopping processes between half of all pairs of 
nearest-neighbour sites can take place. This is not yet full parallel dynamics as desired e.g. 
in traffic-fiow problems, but we will still term it 'parallel'. This model has been considered 
before in the case of deterministic uni-directional motion on the chain [9]. Here we will 
treat the general case where the particles hop with probabilities p and q to the right and 
left, respectively. This contains the deterministic as well as the sequential dynamics of the 
master equation as special cases. The latter is obtained if all probabilities tend towards 
zero. This limiting case corresponds to the Hamiltonian limit in two-dimensional statistical 
physics. In fact, there is a close relation to that area, since the parallel dynamics can be 
formulated as an asymmetric six- vertex model with additional boundary terms. 

It turns out that this general model has properties which are quite similar to those of the 
sequential limit. One physical distinction is a stationary density which alternates from 
site to site and which results from the boundary terms combined with the structure of 
our parallel dynamics. The matrix-product groundstate, which also exists here, has a 
corresponding sublattice structure. The mechanism for the stationary state is the same 
as already encountered in the deterministic limit p = 1, g = [9]. In each time step 
(corresponding to the action of one row of vertices in the vertex model), the sublattices 
exchange their role, so that after two steps the state is reproduced. It should be mentioned 
that simple (scalar) product states in two-dimensional models were studied already some 
time ago. For example, a homogeneous state of that type was found for special cases of 
eight- vertex models with fields [10]. Alternating states were considered in [11, 12] for the 
case of IRF models which are the dual formulation of vertex models. In the context of 
spin systems, such states are called 'disorder solutions' and usually result from competing 
interactions (see [13] and references therein). 

To obtain the stationary state, one may try to find a finite-dimensional representation of the 
four operators defining the algebra and the corresponding vectors satisfying the boundary 
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relations. Here we present a scalar product state (i.e. a one-dimensional representation) 

and a representation in terms of two-dimensional matrices. These representations exist 
for special submanifolds in the parameter space. The general situation, however, is more 
easily handled if one notes that there exist representations where the quantities on one 
sublattice may be eliminated such that one can lift any representation of the algebra in 
the sequential limit with suitably 'renormalized' parameters to a representation of the 
algebra for parallel dynamics. This feature has also been discovered independently in the 
context of a dynamics where the updating is done step-by-step from one end of the chain 
to the other [14] ^). Using this connection, we can take over results and methods for the 
sequential case, e.g. to calculate the current or the density profile in the general case. In 
this way, one obtains a rather complete picture. 

The paper is organized as follows. In Section 2 we introduce the model, discuss the vertex 
formulation, the Hamiltonian limit, and set up relations for the matrix- product state. In 
Section 3 we find one- and two-dimensional matrix representations for the algebra arising 
from the matrix-product ansatz which we use in Section 4 to compute various physical 
quantities. The mapping to the sequential case is described in Section 5 and used in 
Section 6 to compute the current for general values of the parameters. Section 7 contains 
a discussion of the results and remarks on some other aspects as the relation to the model 
in [14] and the integrability of the model. Several appendices contain some supplementary 
material. 

2. Model and matrix-product ansatz 

The diffusion of particles with hard-core repulsion is a stochastic process on a lattice which 
we here choose to be one-dimensional with sites where is even. Each site can have 
two states: It can either be empty or it can be occupied by one particle. Particles can 
hop along the bonds of the lattice onto empty sites, to the right with probability p and 
to the left with probability q. At the left and right boundaries particles can be added or 
extracted. Particles are added to an empty leftmost site with probability a and removed 
from it with probability 7. At the right boundary they are extracted with probability 
P and added with probability 5. All theses processes can take place simultaneously, but 
updates are performed in two steps in order to permit at most one hopping process at 
each site in any time step. Particles are removed and added at the boundaries during the 
first time step. During this time step they can also hop along the bonds connecting the 
even sites with the odd ones to their right. During the second time step they can hop only 
along the other bonds. Because of the restriction that a site may be occupied by at most 
one particle, this is a non-trivial many-body problem. 

For these update rules the time transfer matrix T which describes the time evolution of 
the probability distribution has the structure T = T2T1 where Ti accounts for all processes 
that can take place during the first time step, and T2 for those of the second one. They 



^) In this paper we will use a different terminology than in [14]. We reserve the term 
'sequential' for a situation where in each time step only one local process can take place, 
but at a random position. 
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are given by 
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where the matrices T, C and TZ describe hopping and particle input and output, respec- 
tively. In a suitable basis they are given by 



r 



/I 













l-q 


P 








Q 


1 — p 





Vo 








ij 



1 — a 



a 



7 

1-7 



7^ 



1-5 
6 



1-P 



(2.2) 

This stochastic model can be regarded as a vertex model with the time evolution operator 
T corresponding to the diagonal-to-diagonal transfer matrix [15, 16]. This is shown in Fig. 
1 where the particles sit on the bonds of the lattice and time evolves upwards. The initial 
state before application of T is given by the configuration of the bonds at the lower edge 
of the shaded region, and the final state is given by that at its upper border. 




Fig. 1: Representation of the hopping processes as a vertex model. 
Only a part of the lattice is shown in the vertical direction. The shaded 
region contains all vertices that contribute to the diagonal-to-diagonal 
transfer matrix T. 



The local update operators (2.2) can then be reinterpreted in terms of the Boltzmann 
weights of all possible vertex configurations as shown in Fig. 2 where the presence (absence) 
of a particle on a bond is indicated by an arrow pointing upwards (downwards). This 
vertex model is somewhat more complicated than the one treated in [17]. Usually, one 
considers the symmetric six-vertex model which is invariant under inversion of the arrows 
(i.e. particle- hole symmetry). For p ^ q the bulk vertices in the first line of Fig. 2 do not 
have this symmetry. More importantly, the boundary vertices in the second line do not 
conserve the particle number, because they correspond to particle injection and extraction. 
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Fig. 2: The Boltzmann weights for the vertex model describing hop- 
ping with parallel dynamics. 

In this paper we shall be interested in the stationary state, i.e. a state that is invariant 
under the time evolution operator T. In the language of the vertex model this corresponds 
to the 'groundstate' of the diagonal-to-diagonal transfer matrix. 

Because of the sublattice structure of the transfer matrix (2.1) we make an alternating 
matrix-product ansatz for the stationary state 

(2.3) 

where A, S, A and B are operators in an auxiliary space. The operators A and A describe 
empty places while B and B encode the presence of a particle. \V) and \W) are vectors 
in this auxiliary space which have to be chosen suitably and have to satisfy the condition 
(1^1 F) 7^ in order for (2.3) to be non-zero. 

The mechanism of [9] (and similarly [12]) assumes that Ti as well as T2 exchange the 
operators A, B and A, B with each other. At the boundaries this gives rise to the conditions 

W£(^)=(W|(|), 'j(|)|v) = (b)|v) (2.4) 



and for the interior one has 



B J \B J \ \B J \B 
After inserting the matrices (2.2), this ansatz leads to 

AA = AA, {1 - q)AB + pBA = AB , 

BB = BB, qAB + {l-p)BA = BA, (2.6a) 
{(1 - 5)A + PB}\V) = A\V) , {5A+ (1 - P)B}\V) = B\V) , 

{W\{{l-a)A + ^B} = {W\A, {W\{aA+{l-^)B}= {W\B. (2.6b) 
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It has been argued in [9] that the relations (2.6) define a consistent associative algebra with 
Fock representation. However, so far the problem of finding the groundstatc of T has just 
been reformulated and not yet been solved. In order to make further progress, one needs 
a representation of the algebra defined by the A, B, A and B with suitable additional 
properties. This will be the subject of Sections 3 and 5. 

The corresponding matrix-product state of the sequential limit is well-known (see e.g. 
[3,18,19]). The sequential limit is the limit of small probabilities, or equivalently the 
Hamiltonian limit in the language of vertex models. In order to be more precise set 

X := px (2-7) 

for X = p,q, a, P, 7, S such that one can make an expansion in powers of p. An immediate 
consequence of the definitions is that 

T=l-ph, C = l-phL, n=l-phR, (2.8) 

with matrices h, Hl and ha that are independent of p and describe the local processes 
with rates p, q, a, P, 7 and S. The transfer matrix T now takes the form 

T=l-pH + 0{p^), (2.9) 

where the Hamiltonian H contains only nearest-neighbour interactions. After a similarity 
transformation, this Hamiltonian becomes the [/g(s'u(2))-invariant Hamiltonian of the fer- 
romagnetic XXZ-Heisenberg model [20] with additional boundary terms Hl and hn (see 
e.g. [19]). Since this Hamiltonian does not have any sublattice structure one can make a 
homogenous matrix-product ansatz. We set A — A — E and B = B — D such that the 
operators with and without hat in (2.3) become equal. Then one imposes the following 
relations [3, 18, 19] (see also [8]) at the boundaries 

{W\hL (^^^ = {w\ (^-^^^ , ^h(^)|^) = -("i^)|^), (2-10) 

and the following algebra for the bulk 

With this ansatz the right term in (2.11) at site x cancels the left term at site x -|- 1 if the 
complete Hamiltonian is applied. Finally, one is left with two boundary terms which are 
cancelled by (2.10). 

After inserting the explicit matrices h, Hl and hji one finds that (2.10) and (2.11) are 
equivalent to 

pDE - qED = D + E, 
(PD-'5e')\V)^\V), (2.12) 
{W\ (aE - ^D) = {W\ . 
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This algebra can be used to compute expectation values in the stationary state efficiently 
by means of recurrence relations (see [3, 18, 19]). The operator C := E + D acts like a 
transfer matrix in the spatial direction. Obviously, (VF|C^|y) is the sum of all coefficients 
of the groundstate. In order to obtain a probability distribution one has to divide by 
this factor. Then e.g. the density profile {tx) is given by the expectation values {tx) = 
{W\C''-^DC^-''\V)/{W\C^\V). 

One of the main motivations for the work to be reported below was to find out how the 
mechanism (2.4), (2.5) is related to the known mechanism (2.10), (2.11) in the sequential 
limit. In the next Section we will first see what the relation is in two special cases. 

3. One- and two-dimensional representations 

First we consider a scalar product state. For this representation of the algebra (2.6) the 
operators A, B, A and B are mapped to real numbers a, 6, a and b and the boundary 
vectors can be discarded. (We do not write the representation map explicitly and below 
we will also use the same notation for the algebra and its representation). After fixing the 
normalization suitably, the condition (2.6b) is equivalent to 

a = S^ — ■y — 3 + b — da — a — d + da , 

(3.1) 

a = /?7 — 7 — /3 + Pa , b = S^ — a — S + Sa. 

The bulk relations (2.6a) are satisfied if (1 — q)ab— (1 —p)ba — 0. Inserting the result (3.1) 
one finds that the probabilities p, a, 7 and 5 have to obey the relation 

{I - q){(3-i - - (3 + (3a){5a - a - 5 + (3a) 

= {I - p){5-i - a - 5 + 5a){5-i - -i - (3 + . ^ ' ' 

Note that this condition is symmetric under the exchange of p with q, a with 5 and j3 with 
7. This symmetry is to be expected because parity (reflection at the centre of the chain) 
effects precisely this exchange of the probabilities. A further expansion of (3.2) is not very 
enlightning, but a few special cases may be worth while mentioning. For 7 = 5 = but 
a 7^ ^ /3 the relation (3.2) specializes to 

{l-q){a + l3-al3)=p-q. (3.3) 

This condition agrees with the one that has been found in [14] for the existence of a 
one-dimensional representation with 7 = 5 = 0. Another interesting case is the sequential 
limit. This is obtained by inserting (2.7) into (3.2) and keeping only the leading non-trivial 
order in p. As expected, one obtains a condition that is equivalent to the known result for 
sequential dynamics, namely eq. (78) in [19]. 

We would like to conclude the discussion of the one-dimensional representation with a 
remark on more general mechanisms than (2.4), (2.5). Firstly, one can get rid of the alter- 
nating structure in space by introducing a four-state model that uses block-spin variables 
for the states on any odd site and the even site to its right. Then one can make a ho- 
mogenous product state ansatz. If one further requires the coefficients of this state to be 
independent of N for all even AT > 2, one finds that the relation (3.2) must be satisfied. 
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Even more strongly, one finds that the block-spin product state can be factorized in the 
product state that we have already found above. In this sense, the mechanism (2.4), (2.5) 
gives rise to the most general scalar product state. 

Now we turn to the more interesting case of a two-dimensional matrix-product state. As 
we have already indicated at the end of the previous Section, the matrices 

C:^A + B, C:=A + B, (3.4) 

will play the role of transfer matrices in space. Therefore, it is desirable to find represen- 
tations where C and C have a particularly simple form. First one observes by summing 
over the columns in (2.5) that [C, C] = 0. Since C and C commute they can be simul- 
taneously brought into diagonal -or more generally- Jordan normal form. If C and C 
are both invcrtiblc, an even more stronger statement holds [9], namely that one can find 
an equivalent representation where C — C. On the other hand, one can check that any 
non-trivial two-dimensional representation where C ^ cC with some constant c leads to 
the condition (3.2) and is therefore not interesting. For this reason we will from now on 
only consider representations where 

C = C. (3.5) 

This choice also ensures that half of the boundary equations (2.4) are satisfied automati- 
cally. 

First we consider representations where C = C can be diagonalized. We fix the normaliza- 
tion of the representation by requiring that Ci^i = 1. So, apart from the matrix elements 
e.g. of A and A (and the vectors | V), |VF)), only 6*2,2 is still free. With this choice we first 
solve the bulk equations (2.6a). One solution is given by 



{A{p-q){p+q-l)+p{l-p)){A{p-q)-p) 
iA{p-q)-p+l)pq 

A{q-1) 




A{p-q)-p+l 

rv {q-l)(A{p-q)-p)A 
{Aip-q)-p+l)p 

A 1 

n A{A{p-q){p+q-l)+p{l-py) 
iAip-q)-p+l)p 



(3.6) 



where the constant A remains free. In obtaining (3.6) we have used the freedom to rescale 
the two basis vectors in the representation space independently to set 74.1,2 = 1- 

Now we have to solve the boundary equations (2.6b). Non-trivial solutions for the vectors 
\V) and \W) exist iff 

det (^(1 -a)A + ^B-Aj=0 = det (^(1 - 5)A + PB - A^ . (3.7) 

The solution of (3.7) can be reduced to a discussion of the diagonal elements since all 
matrices in (3.6) are upper triangular. The relevant solution is given by the vanishing of 
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the (1,1) -element in the left matrix and the (2,2)-element of the right one. From this one 
first finds that the constant A is given by 



^ =((1 - p) {pQ'j + q6'y (l-p-q) +pP (p + 'J {I - p - q)))) 

{qS (1 - g) (g (1 - a) + a + 7) + p(3 (1 - p) (p (1 - 7) + a + 7) + 

pg((l -2(5-/3)7+ (1 -2/3-5) a) + pq {{qa + p-f) {/3 + 5 - 1) + q^ 5 + pa (3))~\ 

(3.8) 

and with this result the following relation between the parameters of the model: 
(l-p)(l-g)(pW-gV<^^) 

+ (l-p)pgV («/3(l-g)' + (l-g) iq{a + (3) - {5a + (3^)) - q{{l - d + j - q)) 
-{l-q)p\a(3 (-f5{l-pf + {l-p){p{^ + 5)-{5a + l3-f))-p{{l-a)(3 + a-p)) 
+ pV («7 ((1 - Sf - (1 - f3f) + /35 ((1 - 7)' - (1 - af)) 
+ pV {f3S {1-af+a^ {1 - a) (1-/3)) 

-pV («7 {1-Sf + pS {1-^f-ap (1-7) (1-5)) 

+ pV 5 (1 - a) (1 - /3) - /a /3 (1 - 7) (1 - 5)) = . (3.9) 

This is the condition for the existence of a two-dimensional matrix-product state. Like 
(3.2) this condition is symmetric under the exchange of p with q, a with S and /3 with 7. 

To specify the two-dimensional representation completely one has to determine the bound- 
ary vectors {W\ and \V) that belong to the solution (3.6), (3.8) and (3.9). They are 
straightforwardly obtained by determining the nuUspaces of the matrices in (3.7): 



1^) 



{(3 + d) {Aip-q) + l-p) 
A'^{l-5-P){p-q)-A{p-q+{l-p)5 + P{l + q- 2p)) + /3 (1 - p) 



\W\ = /^(l - {Mp-(1) (pQ - 7P(1 -p) - - q)) + lP^ (l -p)} - ^ocq^ (1 - ^ 
V pq{a + ^)iAip-q) + l-p) J' 

(3.10) 

Inserting herein (3.8) one can check that is neither zero nor singular on the manifold 

given by (3.9). In passing we note that the case of symmetric diffusion in the bulk of 
the chain (q = p) is not permitted because in this case the normalization vanishes, i.e. 
{W\V) = 0. 

In order to illustrate the meaning of (3.9) we are now going to look at a few special cases of 
it. Firstly, consider deterministic hopping exclusively to the right, i.e. p = 1 and q = 0. One 
observes that in this special case the condition (3.9) is satisfied automatically, indicating 
that one can find a two-dimensional representation for arbitrary boundary probabilities a, 
/3, 7, 5 if the bulk rates are p — 1, q — 0. This is consistent with [9] where a two-dimensional 
representation was written down atp=l,g = 7 = 5 = and arbitrary probabilities a 
and /3. However, it should be noted that we may not directly set p = 1 in our solution, 
because then one finds ^ = from (3.8) and singularities appear in (3.6). 



8 



Another interesting case is obtained when ^ = 6 = 0, but p ^ 0, a ^ ^ p. Here (3.9) 
simphfies considerably: 



(1 - g){(l - p - q)a(3 + q{a + (3)} = q{p - q) . (3.11) 

This is a hyperbola in the Q;-/3-plane which intersects the axes at {p — q)/{^ — q)- 

A final case of particular interest is the Hamiltonian limit. This limit of the condition (3.9) 
is obtained by inserting (2.7) and keeping only the leading order in p. Like in the case of 
a one-dimensional representation we recover the known result, namely eq. (81) of [19]. 

One might think that the condition (3.9) arises from the special mechanism (2.4) and 
(2.5). However, for 7 = 5 = we have checked that the solution obtained above is the 
most general two-dimensional one in the following sense: As mentioned at the end of the 
discussion of the scalar product state one eliminates the alternating structure in space by 
introducing a four-state model that uses block-spin variables for two neighbouring sites. 
Then one can make a homogenous matrix-product ansatz with four independent operators 
for each two-site block. One further requires the matrix representation of the operators to 
be independent of the chain length N. Then already a comparison with the groundstate 
of the transfer matrix T at N e {2, 4, 6} leads to the condition (3.11). 

There also exist representations where C = C* is not diagonalizable. We do not discuss 
them explicitly here because we are not going to make use of them in the sequel. The 
interested reader may find some results in appendix A. We mention that non-trivial Jordan 
forms give rise to power laws with positive exponents in the correlation functions (see e.g. 
[8]), in particular for two-dimensional matrices to a density profile that is linear in space. 
Physically, this is related to phase transitions in the system. 

Finally, we turn to the question that was one of our motivations for constructing the above 
matrix-product state, namely how the state given by (3.6), (3.8) - (3.10) behaves in the 
sequential limit. First one observes that the matrices (3.6) have the property 

A-A = B-B= = q. 5)1. (3.12) 

l-p+{p-q)A 

The function g is obtained by inserting (3.8) and too complicated to be presented here. Ob- 
viously, an analogous results holds also for the one-dimensional case because the difference 
of any two numbers trivially is a multiple of the identity. 

Having observed the property (3.12) for the two-dimensional representation, it is natural 
to restrict one's attention to those representations of the algebra (2.6) where A — A and 
B — B are represented by the same multiple of the identity. As we will see soon, such a 
condition enables one to make full use of the machinery developed for the algebra (2.12) 
in [18] and [19]. 

However, let us for the moment concentrate on the one- and two-dimensional representa- 
tions. Using (3.12) to eliminate e.g. A and S, the conditions (2.6) turn into 

{W\ {aA--tB) = g{W\, 

iPB - 5A) \v)=gil-p-5) \V) , (3.13) 
pBA -qAB = g (((1 - q)B + (1 - p)A) . 
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This looks already very similar to (2.12). In order to establish the correspondence man- 
ifestly, one inserts (2.7) into the two-dimensional representation (3.6), (3.8) - (3.10) and 
also into (3.13). Then it is straightforward to check that the limits p — of all four 
operators A, B, A, B exist and are non-zero. Furthermore, one finds that A ^ A 
and S — > S as is to be expected because the alternating sublattice structure must van- 
ish in the sequential limit. The boundary vectors with the normalization as in (3.10) 
have the property \V) = p {\V^^^) + 0{p)) and \W) = {\W^^')) + 0{p)) with non- 
vanishing vectors and Finally, the function g appearing in (3.13) behaves 
as 5f(]?, g; a, /3, 7, 5) = pg{p,q;a, (3,^,6) + Oi^p"^) with a non-vanishing function g of the 
rates in the sequential limit. Now one inserts all this into (3.13) and keeps only the lead- 
ing orders in p. Identifying finally limp^o ^ —'■ and limp^o —'■ one recovers 
precisely (2.12). Thus, we have shown for the one- and two-dimensional representations 
constructed above that the matrix- product mechanism (2.4), (2.5) leads precisely to the 
known mechanism (2.10), (2.11) in the sequential limit. The crucial step to this end was 
the observation (3.12). 

4. Physical quantities from the representations 

In this Section we present results for quantities of physical interest such as the bulk density, 
the current and the correlation length. The computations are based on the one- and two- 
dimensional representation of the previous Section. In order to present the results in a 
compact form it will be convenient to introduce an abbreviation which (as we will explain 
in more detail later) is related to the variables k,± used in [18, 19]: 

^±{x,y) = ^( -x{l-q) + y{l-p)+p-q 

2^ ^ (4.1) 

± ^ (-a;(l -q)+yil-p)+p- qf + ^xy (1 - q) {1 - p)) , 

where either a and 7 or /3 and 5 will be substituted for the arguments x and y. To 
understand the meaning of this quantity consider the special case 7 = 5 = 0. First one 
notes that K+(a;, 0) = —(1 — q) + {p — q)/x and K_(a;,0) = 0. Thus, ioi p ^ q and zero 
boundary probability y, the expression K+(a;,0) is basically the inverse of the probability 
X whereas k_(,x, 0) is trivial. Turning on y > means that one not only injects particles 
e.g. at the left boundary but also extracts them again. Now the variables lij^{x^y) can 
be thought of as effective particle input and output rates. They are non-trivial because 
particles can be e.g. injected at the left boundary, then diffuse some time in the interior 
before they return to the left boundary where they are then extracted again. If this 
interpretation of the K_(_(x,y) is correct, they should be related to the ratio of x and y if 
hopping in the bulk is sufficiently suppressed. Indeed, in the limit p — > 0, g — > one finds 
from (4.1) that K-|-(a;,y) y/x. 

It should be noted that wherever the abbreviations k± are used one has to be careful with 
the range of validity of the result due to the ambiguity of the sign of the square root in 
(4.1). In order to be on safe grounds we restrict to p > q. The results for the case p < q 
can be recovered using parity, i.e. by exchanging p with g, a with 5 and j3 with 7. 

Now we turn to the computation of physical quantities. We shall be interested only in 
the thermodynamic limit A/^ — > 00. In this limit, the computation amounts to taking 
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appropriate matrix elements of the relevant operators as the following argument shows. Let 
Ci,i be the largest of the diagonal matrix elements of C . Then in C'^\V) and in (VF|C"" the 
zth basis vector in the auxiliary space will yield the dominant contribution if n respectively 
m is large. Thus, only the ith matrix elements of the operators corresponding to the 
quantities of interest will contribute in the thermodynamic limit. Below, we will refrain 
from presenting results that are obtained immediately by inserting the matrix elements of 
the one- and two-dimensional representations, but will instead present a formulation using 
the /5± that we have just introduced. These variables are used because we believe that the 
results given in terms of them will be valid also off the manifolds in the phase diagram 
accessible by the one- and two-dimensional representations (see also Section 6). 

We first compute the density profile which is given by {t^) — {W\C^~'^BC^~^\V) / 
{W\C^\V) for X odd and by {%) = {W\C''-'^BC^-''\V) /\w\C^\V) for x even. These 
profiles are independent of x for any scalar product state. Their values can be written 
down immediately for the product state (3.1). The only non-trivial computation is to 
show that the resulting expressions are equivalent to 

= 1 . (fx) = 1 ^I'^i' » (4.2a) 

1 -g + Av+(/3,d) l-p + K+(/3,5) 

and also to 

(Tx) = ^^^^ . , (f.) = r (4.2b) 

l-p + K+(a,7) l-5 + K+(a,7) 

on the manifold given by (3.2). 

For the two-dimensional representation (3.6), (3.8) and (3.10), the density profiles depend 
on the spatial variable x. However, in the bulk of the system, the densities become constant 
and follow easily from the observation made above: (tx) ~ Bi^i/Ci^i and (r^) ~ Bi^i/Ci^i 
for 1 <^ a; <^ A^. For the two-dimensional representation (3.6) and (3.8) two cases have be 
distinguished: 6*2,2 > C*!,! which we denote by T and Ci^i > 6*2,2 which we denote by TP. 
This corresponds to two different regions in the phase diagram, see Fig. 3 below. 

The results obtained from the two-dimensional representation for the bulk (1 -C a; -C N) 
can be identified with the result (4.2) on the manifold given by (3.9). One finds that 
(4.2a) is valid in case I and (4.2b) applies to case II. Note that we recover the result for 
sequential dynamics (eq. (100) of [19]) in the limit of small hopping probabilities. Like in 
[19] we regard the fact that we obtain the same result from the one- and two-dimensional 
representations as an indication that the expressions (4.2) may be valid throughout the 
corresponding phases. For the special case 7 = 5 = they simplify to 

{rx) = l-^^-^, {%) = -^{tx), (4.3a) 
p-q 1- (5 

(T.) = :r^^, (f,) = (l-a)(r.)+a, (4.3b) 
1 — a p — q 

for case I or II, respectively. 

Next we discuss the current. The current can be computed in the interior of the system at 
those places where hopping processes are possible during the next time step and is given 
by expectation values of the operator 

J = pBA-qAB. (4.4) 
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During the first time step it can also be computed at the boundaries where it is given by 
expectation values of the operators Jl = olA — and = f3B — 5 A respectively. 

First, we use this to compute the current from the scalar product state (3.1). The result 
obtained directly from (3.1) can be written on the manifold (3.2) either as 

or as _ 

J _ ^ {p - q) - {a + - q + K_ (a, 7)) 

(p-?)(l-«-7) • ^ 



It is easy to specialize this result to 7 = 5 = after recalling that K_(a;,0) = 0. In this 
special case e.g. (4.5a) turns into J = (3 ^^^z^y^^r^y^ . 

Analogously the current can be computed from the two-dimensional representation. In 
order to obtain the result for the thermodynamic limit one now has to compute it by 
taking the matrix elements Jl = {a{W\A — ^{W\B)- / {{W\C)- at the left boundary, Jr = 

(^PB\V) - SA\V)^ I {C\V))^ at the right boundary and J = [pBA - qAB^ /Cl^ in the 

interior where i is chosen such that Ci^i is maximal. Inserting (3.6), (3.8) and (3.10) here, 
one finds the simplest result at the left boundary for * = 1. The current is equal to (4.5a) 
in case I and to (4.5b) in case II on the manifold where (3.9) is valid. In the limit of 
small probabilities this result goes into the one obtained for sequential dynamics, namely 
eqs. (4.10) and (4.11) of [18]. It is not a coincidence that we have obtained the same 
expressions for the current from the one- and two-dimensional representations, because (as 
we will show in the Section 6) they are valid throughout the corresponding phases. 

Finally, we turn to the correlation length <^ for the two-dimensional representation which 
is related to the eigenvalues of C via exp(l/C) = 6*2, 2/(^1, i- It is given by 

/ 1 \ ^ ^+(a, 7) (1 - P + ^+(/3, 5)){l-q + K+{(3, d)) 
^\CJ ^+(/3,5)(l-p + ^+(a,7))(l-g + ^+(«,7)) ^ ■ ^ 

on the manifold (3.9). To be precise, we have been able to establish the validity of (4.6) 
analytically for either 7 = or 5 = 0. For both rates non-zero we have performed a careful 
numerical verification. As expected, the result of [19] (eq. (102) loc. cit.) is recovered in 
the limit of small probabilities. Whether (4.6) is also valid off the manifold specified by 
(3.9), and if so, in which areas, is beyond the scope of this paper. 

The result (4.6) can also be used to translate the inequalities between 6*2,2 and Ci,i into 
inequalities between K^{a,j) and /5_|_(/3, 5). On the manifold (3.9) one finds numerically 
that k+(q;,7)k+(/?, 5) > (1 — p){l — q) ^). This can be used to infer from (4.6) that 
6*2,2 > C*!,! (case I) iff «+(/?, 5) > k+(q:, 7) and vice versa. 



^) According to [19] one expects that K,^{a,j)K+{(3,d) > (1 — p){l — q) holds for any 
finite-dimensional representation whereof infinitely many should exist. This inequality can 
be obtained by applying the argument of the next Section to the corresponding inequality 
of [19]. 
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5. Mapping onto the Hamiltonian limit 

In this Section we will impose a condition similar to (3.12) for general values of the pa- 
rameters and use this to map the problem for parallel dynamics to the one with sequential 
dynamics. 

After fixing the normalization of the representation of the operators A, B, A and B 
suitably (but different from the one used in Section 3), we will from now on concentrate 
on representations where 

B - B = A- A^l (5.1) 

holds (this condition has also been proposed in [14]). This choice may be a restriction 
in the sense that representations of the algebra (2.6) might exist that are not equivalent 
to one where (5.1) holds. We will now show that (5.1) can be used to lift the Fock 
space representation of the algebra (2.12) to a representation of the algebra (2.6) and 
therefore one can (at least in principle) compute stationary expectation values for arbitrary 
probabilities a, 7, 6, p and q. Even if other representations of the parallel algebra should 
exist, they would have to give rise to the same groundstate that can also be obtained from 
a representation where (5.1) is valid. Thus, from a physical point of view it is completely 
sufficient to study only representations satisfying (5.1). 

We attempt a mapping to the sequential algebra (2.12) by identifying E with a suitable 
linear combination of A and A and D with another suitable linear combination of B and 
13. This leads to the following ansatz for the operators A, A, B and B in terms of E and 
D: 

A = nEE + el, A = ueE - {1 - e)l, 

B = noD - dl, B = nDD + {l-d)l, 

where the free constants n e and n d reflect the freedom of choice of relative normalization 
of the representations and the constants e and d correspond to the points of identification. 
Inserting this ansatz into (2.6a) one obtains precisely one independent relation between E 
and D: 

nEnoipDE - qED) =nEE {d{p - q) + {I - p)) + noD (-e(p - q) + {I - q)) 

(5.3) 

+ ed{p — q) + e(l — p) — d{l — q) . 

In order to be able to identify this with (2.12) the constant term in (5.3) must vanish. 
This is ensured by chosing 

e = (54) 

Also the linear term in E must appear with the same coefficient as the linear term in D 
on the r.h.s. of (5.3). This is the case if 

d{p-q) + {l-p) {d{p - q) + {1 - p)f 

Using (5.4) and (5.5) the relation (5.3) now reads 

n /'f"''\+/'7' 0>PE-,EZ)) = E + P. (5.6) 
{l-p){l-q) 



13 



Inserting the ansatz (5.2) into the boundary equations (2.6b) leads to 



The equations (5.6) and (5.7) are identical to the algebra of the sequential limit (2.12) if 
one identifies the rates in that limit as follows 

_ d(p-q) + (l-p) ^ d(p-q) + (l-p) 

' = ^^ (1-P)(1-.) ^' ^-^^ (1-P)(1-.) ^' 

" = "^ l-ae-7c^ "' ^="^ l + fe + /d-(^ + /3) ^' ('-^^ 

^1 1 , 

7 = ''^^:; ;7 ) o = uez ? ^r, — r? t:^^- 

' 1 - ae - ^d ' 1 + 6e + l3d - {d + P) 

Here, the constants e and no are fixed by (5.4) and (5.5). The constant ue remains 
free and refiects the freedom of normalization of the algebra in the sequential limit. It 
is possible to choose ue such that the 'renormalized' rates in the bulk are equal to the 
hopping probabilities, i.e. p = p and q = q. If one does not fix in a computation of a 
physical quantity, it has to disappear from the final result. Also the constant d is still free 
in (5.8), but for explicit computations it will, in contrast to ue, be fixed to a convenient 
value that ensures ue = no- Then, the operators C for the parallel and sequential case 
may be identified with each other. 

Prom (5.8) one sees that in the sequential limit the renormalized parameters differ from the 
initial probabilities only by a factor he, which also equals n_D because of (5.5). Therefore, 
choosing ue = 1, both sets of parameters and also the matrix-product states become 
identical. 

One can also check that the conditions of Section 3 for having a one- or two-dimensional 
representation can be recovered from known results. All one has to do is insert (5.8) 
where ue and d are kept as free parameters into the conditions for the sequential limit 
[19]. The (i-indepcndent factors of these renormalized conditions are precisely (3.2) and 
(3.9) respectively. So, the simple scalar product state as well as the two-dimensional 
matrix-product state of the general case can be obtained by a simple 'renormalization' 
process from the ones of the sequential limit. The same holds also for higher-dimensional 
representations whose existence was shown in [19]. 



6. Current and phase diagram 



Now we show how to use the mapping of the previous Section to compute quantities of 
physical interest. 

According to [18, 19], the phase diagram in the sequential case is given by a single function 
of the rates at either boundary, i.e. by K^{a,^) and k^{P,S). The same will hold in 
parallel case, if the quantities (5.8) are inserted into the functions k+. This statement has 
already been partially verified in Section 4 where the variables K±{x,y) have been used. 
To be more precise, the parallel case will be descibed by functions K±{a,^) = K±{a,'^) 
and K±(/?, 5) = K±(/3, 5) which are obtained from the k,± of the sequential case [18, 19] by 
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choosing the constant d in the mapping (5.8) such that no = tie- The function k± is given 

by 

where 'K±{x,y) was defined in (4.1). Recall from the beginning of Section 4 that K^{a,'j) 
and K+(/3,5) may be regarded as effective input/output rates at the boundaries. Recall 
also that the k± and k± should be used only for p > g and that the results for the case 
p < q can be obtained by applying parity. 

It is straightforward to compute the currents for the general case with parallel dynamics 
using the results of [18] for the sequential case. The operators needed for that have already 
been given in (4.4) for the bulk and for the boundaries below that equation. Using (2.6b) 
and the ansatz (5.1) one can check that 



_ {W\ {aA-^B)C^-^\V) _ {W\C^-^\V) _ 
^ ~ {W\C^\V) ~ {W\C^\V) 

" {w\c^\v) " ■ 



(6.2) 



This means that the current is simply given by ratios of the normalization constants 
(VF|C^|F) for different N if the condition (5.1) holds. Thus, in order to compute the 
current in the thermodynamic limit, the behaviour of such ratios has to be studied for 
large N. 

One finds from (5.2) that 

C = A + B = nE{E + D) + {e-d). (6.3) 

In [18] the operators D and E have been expressed in terms of the creation and annihilation 
operators F and F''' of a g-deformed harmonic oscillator by setting 

F + 1 + 1 

D = 4^ , E = . (6.4) 

p-q p-q 



Inserting this into (6.3) leads to 



C = nE Z: ^ (6.5) 

p-q 



with 



A = 2 + (e-rf)^=^i=i=i=. (6.6) 
riE 



First we consider the case /?_|_(q;,7) < 1 and k+{/3,5) < 1. According to [18] this is the 
maximal current phase and to obtain J, the same computation as there can be used. In eq. 
(3.39) loc. cit. the coefficient 2 of has to be replaced by A because of the extra constant 
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in (6.5). Tracing the effect of this modification one finds that for parallel dynamics the 
current in the maximal current phase is given by 



j=J_tzl= V^-V^ , (6.7) 

For small p and q the result (4.6) of [18] is recovered. 

This result can be checked [21] using the one-dimcnsional representation because the line 
(3.2) touches the region where (6.7) is valid (see Fig. 3 below). 

The currents for the other two phases can be obtained from similar modifications of the 
computations in [18]. First we consider the case K^{(3,d) > /i_j_(Q;,7) and /i_j_(/3, 5) > 1. 
The relation (4.7) of [18] remains valid with the new parameters and reads 

{W\C^-^F\V)^K+0,6){W\C^-^\V) = K+{/3,6){W\C^-^\V). (6.8) 
One can further rewrite C as follows using (6.3) - (6.6) 




p — q P — Q 

- 5E) + ^ + ^+\^-^> + tip 

5 \ p-q p-Q 



(6.9) 



Making first use of (6.9) and then of (6.8) one finds the following recurrence relation for 
the normalization constants (generalizing eqs. (4.8) and (4.9) of [18]) 

{W\C^\V) 

=y ^:^{(/^+^+?(^-2) -p + q) (W^|C^-'|F)+ (^+5) (WlC'-^FlV)} 
^{[p + S + S{X-2)-p + q) + (P + S) K+{P, S)] {W\C''-'\V) . 

(6.10) 

Inserting this result into (6.2) one recovers after a straightforward computation the result 
(4.5a) that we already obtained from the one- and two-dimensional representations. The 
result for the other phase k_|_(q;,7) > /t_|_(/3,(5) and /t_|_(Q;,7) > 1 can be obtained simply 
by replacing P by a and 5 by 7 and coincides with (4.5b). 

These results can be used to draw the phase diagram which turns out to be essentially 
the same as the one for the sequential limit [18, 19]. It is shown in Fig. 3 in terms of the 
variables k+(x, y), using a logarithmic scale. There are at least three distinct phases which 
are characterized only by the values of K^{a, 7) and S): 

I: K+{(3,d) > K_)_(a,7), Kj^{l3,d) > 1 (high density) 
II: K+(a,7) > Kj^{(3jd), ^+(0,7) > 1 (low density) 
III: /?_|_(a,7) < 1, K^{(3,d) < 1 (maximal current) 

The density in phases I and II is approximately constant in the bulk. On the coexistence 
line k+(q;,7) = ^+{^,5) which separates these two phases the density profile is linear in 
space. 
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Phases I and II can be mapped onto each other using parity and particle-hole symmetry. 
This 'duality' transformation exchanges ^ 1 — f}v+i-a;, a ^ (3, 'j ^ S and keeps p and 
q unchanged. In particular K_|_(a,7) is exchanged with k^{P,S) from which one obtains 
the mapping between the two phases. 

The dotted line in Fig. 3 shows the condition (3.2) for the existence of a onc-dimensional 
representation which can be cast in the form K+{a, 7)k+(/3, 6) = 1. The figure also shows 
the condition (3.9) for having a two-dimensional representation with two choices of p and 
q. In both cases «_(_(«, 7)k_|_(/9, 5) = p/q is verified numerically. For 7 = 5 = 0, it can 
also be shown analytically that this condition is equivalent to (3.9). This suggests that 
the condition for having a two-dimensional representation is a hyperbola in the K+{a, 7)- 
K_i_(/9, S) plane. By changing the value of p/q this hyperbola can be swept over the entire 
region above the dotted line marking the one-dimensional representation. 

102 



10 



K+{P,S) 1 



10-1 



10-2 10-1 1 10 10^ 

k+(q!,7) 

Fig. 3: Phase diagram of the kinetic model. The dotted line shows 
the condition (3.2) for a one-dimensional representation. The condition 
(3.9) for a two-dimensional representation is shown for two choices of 
p and q: The symbol 'o' is for p = 0.75, q = 0.25 and the symbol '•' is 
for p — 0.6, q — 0.4. 

For a more detailed discussion of the case 7 = (5 = we refer to [21]. 

Our results for the current and densities are summarized in Table 1. The correlation length 
is given by (4.6), but probably not throughout the entire phases I and II. Presumably one 
has to choose at least k+(q;, 7)k+(/3, 5) > 1 because on the line ^^{a, 7)k+(/?, 5) = 1 there 
is a scalar product state and the correlation length diverges as one approaches this line. 

All these results have the correct behaviour under the duality transformation which ex- 
changes the phases I and II. 
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Phase 
I 

II 
III 



/3 



a 



J 



(p-q) (l-P-S)^ 

jp-q) - (n + -;)(! - g + /7_(n, 7)) 
(p - (1 - a - 7) 
yi - g - y/l -p 



jr.) 



l-q + K+{f3,S) 
1 - P + «+(«,7) 



jr.) 



l-p + K+{f3, S) 
1 -q + K+{a,^) 



Table 1: Summary of our results for the current J, the bulk density (tx) 
on the odd sublattice and the bulk density {tx) on the even 
sublattice in the three phases. 



We mention that the result for the current can also be obtained by using the bulk-densities 
in the mean-field formula J = p{tx){1 — (tx+i)) — q{l — (Ta;))(''a;+i)- In regions I and II 
this can be shown analytically for 7 = 5 = 0. In region III it can be verified numerically. 



In appendix B we show how to use (3.12) with g = 1 to compute correlation functions for 
the case of symmetric diffusion p — q where the results of this Section and Section 4 may 
not be applied directly (compare also the corresponding remark in Section 3). For this 
special case, the algebra is much simpler. Appendix C describes how one can use the same 
representation to compute physical quantities on finite lattices with system sizes up to a 
few hundred sites. This can be used to check the validity of our results in Table 1 as well 
as the correlation length (4.6) for K^{a,j)K^{l3,S) > 1 numerically. 



7. Discussion 



We have considered the diffusion of hard-core particles between two reservoirs for a par- 
ticular kind of parallel dynamics. The results show that the stationary state has a similar 
matrix-product form as in the sequential case and can actually be obtained from that 
limit. Therefore the physical properties, in particular the phase diagram, are also similar 
although the formulae are more involved. 

In [14] it was found that the same is true for yet another type of dynamics, where the 
stochastic motion takes place step-by-step along the chain. This can be visualized nicely 
in the vertex-model picture as shown in Fig. 4. The processes take place in a diagonal 
strip of the lattice and the time-evolution operator is seen to be the usual row-to-row 
transfer matrix Trow of the vertex model, with a shift in the numbering of the upper row of 
variables and additional boundary vertices. Using the same exchange mechanism at each 
vertex as in Section 2, one finds that the pair of operators A and B only appears in the 
intermediate steps and the matrix-product state becomes a homogenous state involving 
only the operators A and B. An independent treatment of this matrix-product state would 
involve computations similar to those that we have presented here. Alternatively, one can 
also directly use our results for the current and correlation length for the updates as in 
Fig. 4 (the densities are only the same on the even sublattice) [14] . 
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Fig. 4: Representation of a stcp-by-step dynamics starting at the right 
end of the system as a vertex model. 

In order to illustrate the connection to vertex models further, let us also present a graphical 
representation of the matrix-product mechanism used here. Denote the pair of operators 
A, B by a '•' and the pair A and 5 by a 'o'. A dotted line between them indicates that 
their product is to be taken while a termination of the line inside them means that the 
vectors (yV\ and \V) have to be multiplied from the left and right respectively. Then (2.5) 
and (2.4) can be depicted as follows (compare also Fig. 3 of [12]): 




(7.1) 



The first identity is strikingly similar to the Yang-Baxter equation (see e.g. Chapter 2.3 of 
[22]) which for vertex models reads graphically 




(7.2) 



In both cases one can 'pull' a line through the vertex and the corresponding quantities 
exchange places. The boundary Yang-Baxter equations, involving the so-called K-matrices, 
are more complicated. This leads us to the question of integrability of the model. As 
already mentioned, the time evolution operator in the sequential limit is a Heisenberg 
Hamiltonian with boundary fields. It is therefore integrable in the sense that it belongs 
to a whole family of commuting operators [23, 24] . One would expect that his also holds 
for the full vertex model. A proof would have to follow the lines of [23 — 25] . In any case, 
the integrability would be more relevant for the time behaviour of the system than for the 
stationary state which we determined. 
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Occasionally, it has been conjectured that there is a connection between the existence of 
matrix-product states and integrability. However, examples of scalar product states like 
in [10] show that this cannot be true in general. Nor does the construction in [26] ensure 
integrability. It amounts to transforming the bulk relation in (7.1) into a commutator as 
in (7.2) by going over from the eigenvector |$) to the operator P = |$)($|. However, this 
does not help in finding other eigenstates because P is a simple projector whose unique 
eigenvector with eigenvalue 1 is |$) while a highly degenerate eigenvalue accounts for all 
other vectors. 

Finally, one may ask if the parallel dynamics considered in this work might be used in 
other related problems. Here, two situations come to mind. One would be a model with 
coagulation and decoagulation in addition to the hopping processes. With a proper tun- 
ing of the parameters, this is a free-fermion problem in the sequential limit [27] and the 
stationary state has been shown to have a matrix-product form with four-dimensional ma- 
trices [8] . As in the present work, there are four matrices involved, two for the homogenous 
matrix-product state and two in the generalization of the cancellation mechanism (2.10) 
and (2.11). Indeed, if the mechanism (2.4) and (2.5) should be applicable to more general 
situations than discussed here, one would in general not expect the differences A — A and 
S — S to be proportional to the indentity and thus one would find the more general matrix- 
product mechanism of [8] in the sequential limit. However, the choice C = C is always 
possible independent of the details of the dynamics. This would lead to a linear relation 
between the matrices in the sequential limit, but one can check that the four matrices 
used in [8] are linearly independent. Thus, whether the problem including coagulation and 
decoagulation admits a matrix-product state also for parallel dynamics, and if so, with 
which mechanism, remains to be investigated. 

The other problem is hopping on a ring with a defect where the rates are modified. For- 
mally, this case is obtained by replacing the product TZ® C oi the boundary matrices by 
a hopping matrix T . This model has already been solved by a modified Bethe ansatz for 
the case of unidirectional deterministic motion everywhere except at the defect [28]. For 
a certain fixed particle density, the stationary state is also expressible in the form of a 
two-dimensional matrix product [29]. One may speculate that our representation (3.6) of 
the bulk algebra may help to solve the general model with a defect. In particular, the 
observation that the differences of the matrices with and without hat are proportional to 
the identity is a pure bulk property and thus one may hope that techniques like those we 
have used e.g. in Sections 5 and 6 may be useful also for systems with a defect. 

Acknowledgments 

We thank H. Hinrichsen, M. Karowski, P. Pearce, N. Rajewsky, V. Rittenberg and A. 
Zapletal for useful discussions. A.H. would like to thank the Deutsche Forschungsgemein- 
schaft for financial support. 



20 



Appendix A. Representation with Jordan form 

In this appendix we briefly discuss a two-dimensional representation where C = C has 
a non-trivial Jordan form. This case occurs if 6*2,2 = 1 and simultaneously = 

because then the normalization (VF|C''^|y) of the groundstate vanishes and the representa- 
tion with diagonal C may not be used any more. We have not been able to fully determine 
when this happens in (3.6), (3.8) - (3.10), but one can check that it is the case if P — a 
and 7 = 5. Using these values in (3.9), one then finds that 

{pa^ - qS^) (1 -p)(l-q)+pq(^(l- p) (1 - 5f - {1 - q) {1 - af) = . (A.l) 

One can choose 

C = e=(j J) (A.2) 

by suitably fixing the normalization of C and C and the relative normalization of the two 
basis vectors in the auxiliary space. The further strategy to determine the representation 
matrices for A und A is now slightly different from the procedure described in Section 3. 
First one solves the boundary conditions (2.6b) as described below (3.7). With the result 
of this computation one can then solve the equations (2.6a). This leads to: 

^_^fqS — S + pS— p 



B \ qa + q6-a-q 

/ ^ ^ \ (A.3) 

■^_l[pd — d+pa— p 



B \ qa — a + pa — q 

with 

B^{a + 6){p + q-l)-{p + q). (A.4) 

Again, some freedom was used to set Ai^2 = in (A. 3). Instead of (3.10) one now has to 
choose the boundary vectors as follows: 

\^^ = {{a + S){{l-S)il-p)-il-a)il-q)}) ' 
^^^^l^-{a + S){il-S){l-p)-{l-a)il-q)r 

One can check that is neither zero nor singular on the manifold given by (A.l). 

A property similar to (3.12) also holds here and has in fact been assumed to simplify the 
computations. 

We should mention that (A.l) is not recovered after inserting the renormalized rates (5.8) 
into its sequential limit. However, it is argued in Section 6 that the coexistence line has a 
more general form than just /3 — a, 'j = 5. Thus, the fact that (A.l) does not reproduce 
itself under renormalization is a hint that one should be able to find two-dimensional 
representations where C is given by (A.2) under more general conditions than (A.l). 
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Appendix B. Symmetric difFusion in the bulk 

Here we discuss the special case of symmetric difFusion in the bulk, i.e. p = q. In this case, 
which was also considered in [15], it is convenient to work with the operator B and the 
spatial transfer-matrix C. Inserting the choice (5.1) into (2.6) leads to 

[B,C] = i^C, 

^ ^ P (B.l) 
{W\ {aC - (a + 7)S) = {W\ , ((/? + S)B - SC) \V) = {1 - (3 - S) \V) . 

Using these relations one immediately finds by either commuting the B to the right or the 
left boundary 



■^(- - 1) - ^^i^) + ^^^<'^ic"i''> • 

(B.2) 

From these two expressions the term (VF|C^ ^W) can be eliminated yielding the density 
profile for the odd sites 

a (/3 + 5) (iV - 1) - (a - 7 5) (x - 1) + ^ {a (1 - /3) + 5 (1 - a)} 

^"""^ " (« + 7) (/3 + 5)(iV-l) + ^{a (l-/3)+5 (l-a)+7(l-5)+/3(l-7)} ' 

(B.3) 

A similar computation with B = B -\- \ yields the profile for the even sites. The ex- 
pression differs from (B.3) only in the last curly bracket in the numerator which becomes 
{a (1 — 5) + 5 (1 — 7)}. In the thermodynamic limit N ^ 00 this difference between (r^) 
and (fj;) disappears and, more strikingly, the densities become independent of the bulk 
probability p. The result after the limit is identical to the one obtained in [30] for the case 
of sequential updates. 

For finite N, the density profile (B.3) is linear in the spatial variable x (as one can 
already see from (B.2)). For p = q the current is related to the density profile by 
J = p ((fa;) — {tx+i)). Inserting the result (B.3) one finds that 

'-''{l-p)ia + ,)if3 + sf ^^-'^ 

for large. Physically, this means that the pumping effect at the ends is not sufficient 
to drive a current through an infinite system. The results (B.3) and (B.4) show that only 
one phase exists which corresponds to the coexistence line in Fig. 3. 

It is straightforward to apply the method used above for computing the density profile also 
to higher correlation functions. 

Although for the case discussed here the groundstate is comparably simple, this is still 
interesting because for p — q our local bulk-transfer-matrix T becomes identical to the 
i?-matrix R of the symmetric six- vertex model. One can parametrize p = u/{u + rj) where 
u is the spectral parameter and 77 a constant. With this parametrization the Yang-Baxter 
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equation (7.2) with spectral parameter u is satisfied and can be used to write down a family 
of commuting stochastic row-to-row transfer matrices for periodic boundary conditions 
(see e.g. Chapter 2.3 of [22]). Thus, the hopping probability p is essentially the spectral 
parameter of an integrable model. One consequence probably is the disappearance of p 
from the density profile (B.3) in the thermodynamic limit. 

For p q, the T-matrices (2.2) can be identified with the ^-matrix of the general asym- 
metric six-vertex model with a suitable spectral parameter u: 



R 



/sm{u + r]) ^ \ 

z'^ sm{r]) X sm{u) 

sin('u) 2;~" sin(r7) 

V sm{u + r])J 



(B.5) 



where ry, x and z are some arbitrary constants. One can show that this also satisfies the 
Yang-Baxter equation (7.2). However, the impact of the boundary vertices in Fig. 2 on 
the integrability remains to be fully clarified. 

Appendix C. Computations on finite chains 

Wc show here how to compute correlation functions efficiently using a representation of 
the algebra where (5.1) holds. Then A and B satisfy (3.13) with g — 1. A convenient basis 
for the Fock space is given by B^\V) and one then has to compute AB^\V). In order to 
describe how this can be done e.g. on a computer let 

\x) :^ B^'IV) , \x,y) :=B''ABy\V). (C.l) 

The bulk relation in (3.13) (with ^ = 1) can be used to commute an operator A one place 
to the right which in terms of the above vectors yields the relation 

\x, y) = ^ {p\x + l,y - 1) - {1 - p)\x,y - 1) - {1 - q)\x + y)) (C.2) 

for y > 0. As soon as the operator A hits the right boundary one can use the second 
boundary equation in (3.13) with ^ = 1 to replace Ahy B: 

\x, 0) = i (/?|a; + 1) - (1 - (/? + 5))\x)) . (C.3) 

By definition, the algebra acts on the states (C.l) as 

A\x) = \Q,x), B\x) = \x+l). (C.4) 

The rules (C.2) - (C.4) are sufficient to express any state ABAA. . .\V) in terms of the 
\x) and one does not need the original algebra anymore. In order to be able to compute 
words (i.e. scalar products with (W\) a final constant Sx := (VF|a;) is needed. It can be 
computed by creating an A at the left boundary using the first boundary relation in (3.13) 
(with = 1). This leads to 

Sx+i ^-{a{W\Q,x)-Sx) (C.5) 
7 
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which amounts to a recurrence relation for the after moving the A to the right boundary 
using the previous relations. Setting sq = 1, one is now able to compute the value of any 
word {W\ABAA ...\V) exclusively from the rules (C.2) - (C.5). 

It is possible to solve these recurrence relations in closed form following the lines of [18]. 
However, also formulae of the type as in [18] are best evaluated numerically using a recur- 
sive procedure. The recipe presented above is sufficient to compute correlation functions 
numerically on finite chains where a length = 100 is no major problem. However, in 
order not to do a computation twice, one should store the expansion of the vectors \x, y) 
in terms of the basis vectors |x). One also needs to be careful with the numerical range 
because some of the numbers grow exponentially with N, e.g. (VF|C^|y) ~ and 
|J-i| >1 can become quite large ). 

Fig. 5 shows a density profile obtained in this manner on a finite lattice with = 200. The 
parameters yield K+{a,j) = 0.912 and K+{f3,d) = 0.801 which to corresponds to a point 
at the top right corner of phase III in Fig. 3. The distinction between the two sublattices 
is clearly visible, so is the infiuence of the boundaries. As a byproduct in this computation 
one finds the current with these parameters for = 200 using (6.2): J = 0.2690. This is 
to be compared to the result J = 0.2679 obtained from (6.7) in the thermodynamic limit. 
Computations like this one provide room for more detailed investigations. 
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Fig. 5: A density profile in the maximal current phase with = 200 
sites and p = 0.75, q = 0.25, a = 0.5, f3 = 0.6, 7 = 0.1 and 5 = 0.2. 
The full line shows the density {t2x+i) and the dashed line the density 




^) An implementation in C which takes care of all these details is available on the WWW 
under URL |http: / /www. physik.fu-berlin.de/~ag-peschel/software/mp. html] . This pro- 
gram does not only compute he current, the density profile and the two-point function on 
a finite chain, but also implements our results in Table 1 as well as the correlation length 
(4.6) and thus provides a simple way for checking their validity numerically. 
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